Application of realistic effective interactions to the structure of the Zr isotopes 
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We calculate the low-lying spectra of the zirconium isotopes [Z = 40) with neutron numbers 
from A'^ = 52 to A'^ = 60 using the lpi/20(?9/2 proton and 2slci0p7/20/iii/2 neutron sub-shells to 
define the model space. Effective proton-proton, neutron-neutron and proton-neutron interactions 
have been derived using **Sr as closed core and employing perturbative many-body techniques. The 
starting point is the nucleon-nucleon potential derived from modern meson exchange models. The 
comprehensive shell-model calculation performed in this work provides a qualitative reproduction 
of essential properties such as the sub-shell closures in ^''Zr and ^®Zr. 
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O ' I. INTRODUCTION 

o . 

The Zr isotopes undergo a clear and smooth shape transition with increasing neutron number. The isotopes which 
^ ^ are displayed in Figure 0(a) span from pure spherical nuclei that can be described in terms of simple shell-model 
configurations to the strongly deformed nucleus ^^■^Zr. Evidence for coexisting shapes has been reported around ^°'^Zr 
In the intermediate region, both ^^Zr and ^^Zr present evidence for sub-shell closure. These isotopes have a 
C remarkably large gap between the ground state and the 2f level, 1.751 MeV in ^^Zr and 1.223 MeV in ^^Zr. They also 

have a relatively low level density below 3 MeV. Empirically, the Zr isotopes are fairly well established. Lhersonneau 
^ ' and collaborators have recently performed careful experimental studies of the Zr isotopes and neighbouring nuclei, 
, and they have made major contributions to the identification of levels in ^^Zr and ^^Zr 

' For comparison the empirical Sr spectra are sketched in Figure 0(b). The Sr isotopes differ from the Zr isotopes 
\^ by only two protons but are qualitatively quite different. Compared to Zr the Sr isotopes have a much smoother 
' behaviour with a stable 0+ — 2+ spacing similar to that observed in tin. Large gaps due to sub-shell closure as 
] observed in ^^Zr and ^^Zr do not occur in Sr. However, as in Zr there is a clear transition from spherical to deformed 
■ shape around N = 60. The low-lying spectrum of ^"''Sr is a nearly perfect rotational band. It is a theoretical challenge 
to describe a sequence of isotopes with such big changes in the structure from one nucleus to another as in Zr. For a 
proper description of the Zr isotopes one has to allow for proton excitations. In particular, the protons seem to play 
G ' a dominant role in the state. Thus, the common choice of inert core has been *^Sr, though with large variations in 
I , the size of model space, truncation scheme and effective interactions, Refs. |^-^. The early calculation by Auerbach 
^ ' and Talmi |^ was carried out with valence protons filling the (1^1/2,0(79/2) oribitals and valence neutrons filling the 
^ ' ( 1*^5/2 ) orbital. In particular if one wants to describe more neutron rich Zr isotopes a larger model space is required 
^ ' as the interaction between (0179/2) protons and (2c?3/2) and (0*77/2) neutrons becomes increasingly important Q. To 
J> ' our knowledge the present work is the first shell- model calculation of the Zr-isotopes which includes nuclei up to ^^Zr 
within a non-truncated p : (1^1/2,0(79/2) n : (1^5/2, 2si/2, 1^3/2,0^7/2,0^11/2) model space and with a fully realistic 
effective interaction. 

^ . We will in the present work perform a systematic shell-model study of the Zr isotopes from N=50 to N=58, and 
will in particular pay attention to the nuclei around the closure of the neutron (1*^5/2) and (2si/2) sub-shells, ^^~^^Zr. 
In several works we have performed thorough analyses of the effective two-body interaction. We have derived shell- 
model effective interactions based on meson exchange models for the free nucleon-nucleon interaction, using many-body 
perturbation theory as described below. The systems that have been studied are reaching from the oxygen region to 
the tin isotopes and the N=82 isotones, Refs. [p0|-Jl3[. For the lighter systems, such as the sd- and p/-shell nuclei, 
we obtained markedly better results for nuclei with one kind of valence nucleons than for nuclei with both kinds. For 
the heavier systems we have so far restricted ourselves to nuclei with like valence nucleons, such as the Sn isotopes 
and the N — 82 isotones. This way we have managed to keep the dimensionality of the eigenvalue problem within 
tractable limits. Further, we have seen the need for establishing confidence in the T = 1 interaction before considering 
systems with both valence protons and neutrons, where the proton-neutron interaction may play a crucial role. In 
fact, the Zr isotopes represent a challenge on both these accounts. 

Let us in terms of a simplified model like the weak coupling scheme, in which the proton-neutron interaction is 
assumed to be weak, see if one can gain insight into the qualitative properties of this interaction. In Figure we 
demonstrate the validity of this scheme by seeing how well it describes properties of 92,94,962j.^ column one the 
empirical spectrum of ^"Zr (dashed lines), which represents the proton degrees of freedom, is put on top of the ^"^Sr 
spectrum (solid lines), which represents the neutron degrees of freedom. This would represent the ^^Zr spectrum in 
the weak coupling limit and should be compared with the empirical ^^Zr spectrum in column two. Similarly the weak 
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coupling spectra ^^Sr+^°Zr and ^ Sr+^°Zr are compared with the empirical ^ Zr and ^^Zr spectra, respectively. Most 
states in ^^Zr and ^''Zr are well reproduced by the weak coupling scheme. The model does however collapse in ®^Zr, 
due to the presumed closure of the 10(5/2 sub-shell which does not have a counterpart in ^''Sr. The fact that the weak 
coupling scheme is fairly successful in describing ^^Zr and ^"^Zr leads us to believe that the proton- neutron part of the 
effective interaction is either rather weak or state-independent. 

The paper is organized as follows. In Sect. 2 we give a summary of the calculation of the effective interaction and 
the shell model. Then the results are presented and discussed in Sect. 3 and conclusions are drawn in Sect. 4. 



II. THEORETICAL FRAMEWORK 



The aim of microscopic nuclear structure calculations is to derive various properties of finite nuclei from the 
underlying hamiltonian describing the interaction between nucleons. When dealing with nuclei, such as the Zr isotopes 
with A = 90 — 100, the full dimensionality of the many-body Schrodinger equation 

H^,il,...,A)^E,^>,il,...,A), (1) 

becomes intractable and one has to seek viable approximations to Eq. (Q). In Eq. (|^), Ei and "ifi are the eigenvalues 
and eigenfunctions for a state i in the Hilbert space. 

One is normally only interested in solving Eq. (|^) for certain low-lying states. It is then customary to divide the 
Hilbert space into a model space defined by the operator P and an excluded space defined by a projection operator 
Q = l-P 

d 00 

p = ^|^,)(V'.| Q= E (2) 

i=l i=d+l 

with d being the size of the model space and such that PQ — 0. The assumption is that the low- lying states can be 
fairly well reproduced by configurations consisting of a few particles occupying physically selected orbitals, defining 
the model space. In the present work, the model space to be used both in the shell-model calculation and in the 
derivation of the effective interaction is given by the proton orbitals lpi/2 and 0.99/2 and the neutron orbitals 2si/2, 
lc^5/2, lc?3/2, O57/2 and O/111/2. 

Eq. (]l|) can then be rewritten as a secular equation 

PH,sP^, = P{Ho + Feff = E,P^„ (3) 

where iJ^ff now is an effective hamiltonian acting solely within the chosen model space. The term Hq is the unperturbed 
hamiltonian while the effective interaction is given by 

00 



with V^g\ Ks"*' Kff^--- being effective one-body, two-body, three-body interactions etc. It is also customary in nuclear 
shell-model calculations to add the one-body effective interaction V^g-' to the unperturbed part of the hamiltonian so 
that 

iJcff = Ho + Kff + Kff + ■ • ■ , (5) 

where Hq = Hq + V^g' . This allows us to replace the eigenvalues of Hq by the empirical single-particle energies for 
the nucleon orbitals of our model space, or valence space. Thus, the remaining quantity to calculate is the two- or 
more-body effective interaction X]t:2 ^cS ■ ^^^^ work we will restrict our attention to the derivation of an effective 
two-body interaction 

VoB = v}^\ (6) 
using the many-body methods discussed in Ref. H] and briefly reviewed below. 
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A. Effective interaction 



Our procedure for obtaining an effective interaction for the Zr isotopes starts with a free nucleon-nucleon interaction 
V^^"^ which is appropriate for nuclear physics at fow and intermediate energies. At present, there are several potentials 
available. The most recent versions of Machleidt and co-workers the Nimjegen group and the Argonne group 
I prt have all a per datum close to 1. In this work we will thus choose to work with the charge-dependent version of 
the Bonn potential models, see Ref. ||l5|. The potential model of Ref. jl^ is an extension of the one-boson-exchange 
models of the Bonn group jl^ , where mesons like tt, p, 77, S, lu and the fictitious a meson are included. In the 
charge-dependent version of Ref. the first five mesons have the same set of parameters for all partial waves, 
whereas the parameters of the a meson are allowed to vary. 

The next step in our perturbative many-body scheme is to handle the fact that the strong repulsive core of the 
nucleon-nucleon potential V is unsuitable for perturbative approaches. This problem is overcome by introducing the 
reaction matrix G given by the solution of the Bethe-Goldstone equation 

G = V + V %-G, (7) 

LU - QTQ 

where co is the unperturbed energy of the interacting nucleons, and Hq is the unperturbed hamiltonian. The operator 
Q, commonly referred to as the Pauli operator, is a projection operator which prevents the interacting nucleons from 
scattering into states occupied by other nucleons. In this work we solve the Bethe-Goldstone equation for five starting 
energies uj, by way of the so-called double-partitioning scheme discussed in e.g., Ref. Q. The G-matrix is the sum 
over all ladder type of diagrams. This sum is meant to renormalize the repulsive short-range part of the interaction. 
The physical interpretation is that the particles must interact with each other an infinite number of times in order to 
produce a finite interaction. 

Finally, we briefly sketch how to calculate an effective two-body interaction for the chosen model space in terms 
of the G-matrix. Since the G-matrix represents just the summmation to all orders of ladder diagrams with particle- 
particle intermediate states, there are obviously other terms which need to be included in an effective interaction. 
Long-range effects represented by core-polarization terms are also needed. The first step then is to define the so-called 
Q-box given by 

pqp = pgp + p(g '^^ G + G ^^ G ^^ g + ..\p. (8) 

\ LU — Hq lu — Hq LU — tlQ J 

The Q-box is made up of non-folded diagrams which are irreducible and valence linked. A diagram is said to be 
irreducible if between each pair of vertices there is at least one hole state or a particle state outside the model space. 
In a valence-linked diagram the interactions are linked (via fermion lines) to at least one valence line. Note that a 
valence-linked diagram can be either connected (consisting of a single piece) or disconnected. In the final expansion 
including folded diagrams as well, the disconnected diagrams are found to cancel out |p^ . This corresponds to the 
cancellation of unlinked diagrams of the Goldstone expansion These definitions are discussed in Rcfs. p^ , p^ . 

We can then obtain an effective interaction Hcs — Hq + ^cff^ ^-"^ terms of the Q-box with 

v;S'(»)^<j+i:i£?{'^s'(..-i)r. <») 

m— 1 

where (n) and {n — 1) refer to the effective interaction after n and n — 1 iterations. The zeroth iteration is represented 

(2) 

by just the Q-box. Observe also that the effective interaction V^g (n) is evaluated at a given model space energy lu, as 
is the case for the G-matrix as well. Here we choose lu = —20 MeV. Moreover, although Q and its derivatives contain 
disconnected diagrams, such diagrams cancel exactly in each order [l9| ], thus yielding a fully connected expansion in 
Eq. (^). Less than 10 iterations were needed in order to obtain a numerically stable result. All non- folded diagrams 
through third order in the interaction G are included. For further details, see Ref. [Q. 



B. Shell model 



The effective two-particle interaction can in turn be used in shell- model calculations. Our approach in solving 
the eigenvalue problem is the Lanczos algorithm, which is an iterative method that gives the solutions of the lowest 
eigenstates. The technique is described in detail in Refs. p^,^. The shell-model code developed by us is designed 



3 



for an m-scheme Slater determinant (SD) representation. Even with a rather restricted single-particle basis the size 
of the shell-model problem grows rapidly with increasing number of active valence particles. Table Q shows how the 
number of configurations in an SD-basis grows with the number of valence particles acting within the (lpi/2 0179/2) 
proton shell and the (2sl(i0(77/20/iii/2) neutron shell relative to the ®*Sr-core. Note that the ^^Zr system consists of 
more than 26.000.000 basis states. Only a few years ago shell-model calculations on such systems were not tractable. 
Even with today's fast computers and effective algorithms this kind of calculation is still rather time consuming. 

The single- neutron energies are taken to be those deduced from ®^Sr in Refs. |2 l,P3|]. I n the literature the single- 
proton energy splitting e(lpi/2) — £(039/2) varies from 0.839 MeV to 1.0 MeV |l7|,^| j24| - |2q| . As the final results show 
little sensitivity to variations within this energy interval we let the Oc/g/-?, single-proton energy relative to lpi/2 be 0.9 
MeV. The single-particle energies used in this work are listed in Table 



III. RESULTS AND DISCUSSIONS 



In this section we present the calculations of the Zr isotopes. Firstly, we analyze some systematics of the even Zr 
isotopes. The results are displayed in Figures I and I (for some selected states) and in mo re de tail in Ta bles ^ and 
XI. Then, we proceed with the discussion of the odd Zr isotopes, in Figure I and Tables ^ and ^Il|. A major aim 
is to investigate the effective interaction that has been derived for this mass region. Since the odd nuclei are generally 
more sensitive to the underlying assumptions made, this may give an even more severe test of the interaction and 
the foundations on which our model is based. At the end of this section we discuss problems concerning the binding 
energies. In order to study the effect of the proton degrees of freedom we have also performed a more restricted 
calculation, considering only valence neutrons with respect to a ^"Zr-core. 



A. Even isotopes 

The experimental spectra of the two nuclei, ^^Zr and ^"^Zr, show very similar features. The shell-model calculation 
does also provide 0^, 2^ and 4^ levels in ^^Zr and ^''Zr with quite similar features, although the three levels are 
too compressed compared to their experimental counterparts. Comparison of the results obtained with ^'^Zr- and 
®^Sr-cores indicates that these levels are little affected by proton excitations. In ^^Zr the calculated 0^, 2^ and 4^ 
levels are still too compressed, though not as pronounced as in ^^Zr and ^'*Zr, while in ^^Zr the calculated spectrum 
is more open than the experimental one. 

The energy of the empirical 3j~ level. Figure |l|(a), is monotonously reduced with increasing neutron number. In ^^Zr 
and ^''Zr the calculated 3j~ level is obtained at about 2 MeV, while in ^^Zr and ^®Zr the 3j~ level is obtained much too 
high at about 3.5 MeV. Due to the extension of the model space from a ^°Zr-core to a ®®Sr-core the 3j~ level in ^^Zr 
and ^^Zr undergoes a considerable lowering, yielding results close to the experimental values. From the occupation 
numbers in Table III it is clear that the 3^ state undergoes a structural change from ^"^Zr to ^^Zr. The 057/2 and 2si/ 
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orbitals start to play a more important role. For comparison, the observed 3j~ levels in Sr, Figure |l|(b), are all situated 
around 2 MeV. Their calculated counterparts arc located too high, at about 3 MeV excitation energy. The structure 
of the 3j~ level in Zr is totally different from the structure in Sr, Figure ^ Let us look in detail into these states by 
comparing the 3i levels in ^^Sr and ^"^Zr. Both nuclei have four valence neutrons. In ^^Sr the dominant configuration 
is [(lrf5/2)"^0ft.ii/2] j=3- whereas in ^'^Zr the dominant configuration is [(lpi/20(79/2)j„=4-,5- (1^5/2) 7 =2+].-'=3"- 
difference can be ascribed to the single-particle energies. In Zr the 3^ state is created by exciting a proton into the 
0(79/2 orbital instead of exciting a neutron into the 0/iii/2 orbital. Because the Oft.11/2 orbital is located very high in 
the single-particle spectrum it is more favourable to excite a proton rather than a neutron. 

The situation for the 5^ level is more stable throughout the whole sequence of Zr isotopes, with empirical values 
between 2.5 and 3.0 MeV. The ^°Zr-core calculations provide level energies that are 1.0 - 1.5 MeV too high, while the 
®*Sr-core calculations give energies too low by about 1 MeV. This state consists predominantly of configurations with 
a proton excited into the 0179/2 orbital and the neutrons remaining in the lowest possible single-particle orbitals. 

As pointed out in the introduction, there are strong variations in the structure from one nucleus to another, refiected 
in for instance the 0+ — 2+ spacing. Qualitatively we reproduce the variation in the 0+ — 2+ spacing quite well, as 
shown in Figure |^, although in ^0-96-^^. ^j^g g^^p jg 200 — 400 keV less than the experimental spacing. 

In spite of clear differences in the experimental Zr and Sr spectra, the shell-model calculation provides rather similar 
results. The calculated Zr spectra are in far better agreement with experimental data than the calculated Sr spectra, 
which may indicate that the core is in a different condition in the two systems. In Sr the core seems to be relatively 
soft, whereas in Zr the two additional protons tend to stabilize the core. 
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1. Proton configurations 



The proton degrees of freedom are crucial in order to describe certain energy levels. For example, the first excited 
state in ^^Zr has strong components of proton excitations. From the occupation numbers in Table IV we see that 
the character of 0^ state is totally different from the 0^ state. The proton parts of their wave functions are almost 
orthogonal to each others. In the ground state the protons are most likely to be found in the lpi/2-orbital, while for 
the excited 0+-state it is more probable to find the protons in the 0g9/2-orbital. However, in ^^Zr the two shell-model 
0^ states have nearly the same proton structure, almost pure {1p\i2)tt configuration. The 0^ levels in ^^Zr show 
similar structure as in ^^Zr, but the state has slightly stronger 059/2 mixing than in ^^Zr. The third 0'*' state in 
both ^®Zr and ^*Zr have (0^9/2) as the predominant proton configuration. The change in the proton configuration 
with increasing neutron number was observed in pick-up experiments by Saha e< al. p7|. 



2. Sub-shell closure 



Clear signs of sub-shell closure are seen in ^^Zr and also in ^*Zr, due to filling of the 1^5/2 and the 2si/2 orbitals, 
respectively. From ^''Zr to ^^Zr the gap between the ground state and the 2^ state is doubled. The experimental 
0^ 2^ gap increases, from 0.919 MeV in ^''Zr to 1.751 MeV in ^^Zr, and the calculated gap is also more than doubled, 
from 0.520 MeV to 1.426 MeV. In '^^Zr the calculated spacing is larger than the experimental one. The experimental 
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0+ — 2+ spacing is 1.223 MeV, and the corresponding calculated spacing 1.463 MeV. The ^^Zr ground state Id^/ 
occupation number is 5.66, and in ^*Zr the lc?5/2 and 2si/2 occupation numbers are 5.76 and 1.87, respectively. 

The total impression of the ^^Zr shell- model results displayed in Figure ^ is disappointing. Only the 2^ level is 
reasonably reproduced. As an alternative to the extremely time and space consuming calculation presented in Figure 
we may close the lc?5/2 orbital and perform a ^^Sr-core shell- model calculation of the system. The results, shown 



in Table XI, are much improved. 



3. E2 transition rates 

Experimental and calculated E2 transition rates are tabulated in Table 0. In order to bring the theoretical results 
into agreement with the measured 2^ — > 0^ transition rates we have employed effective charges of 1.8e and 1.5e 
for the protons and neutrons, respectively. These values are consistent with the effective charges obtained in Refs. 
P, p6|p^ , however a bit overestimated compared to the fitting to data on ^^Mo and ^^Zr done by Halse in Ref. To 
be mentioned, the calculation of effective charges based on perturbative many-body methods p9| ] gives much smaller 
values, l.le and in the range 0.5e — 0.7e for the proton and neutron effective charge, respectively. 

We adopt Halse's effective value of the oscillator parameter h = 2.25 fm. The value was choosen by reference to 
measurements for the radii of the single-particle orbitals in ^^Sr and of the charge distributions in ^^~^^Mo, Refs. |2|] 
and 0. 

With effective charges and the oscillator parameter as described above, the transition rates between yrast states are 
fairly well reproduced. In the former discussion we have focused on the 0^ state, in particlular its proton structure. 
In ^^Zr and ^^Zr we totally fail in reproducing the transition rates involving the Oj state. The experimental transition 
rates between O2 and 2^ in ^^Zr and ^'^Zr are relatively strong, 14.3(5) and 9.3(4) W.u., respectively, whereas the 
calculated transition rates are two to three orders of magnitude smaller. Similarly, in ^^Zr there is an experimental 
transition rate between Oj^ and 2^ with strength 51(5) W.u. The corresponding calculated transition rate is negligible. 



In fact, the occupation numbers in Table IV show that the 0;^ state in ^^Zr has similar proton structure to the Oj 
state in ^^Zr and ^"'Zr. Leaving out the neutron contributions, as done in column 6 of Table by setting the effective 
neutron charge equal to zero, we observe that the proton part of the wave function contributes more to the transitions 
involving the excited 0"^ states than what it does to the other transition rates. The contribution is however far 
from sufficient and there is a cancellation effect between the proton and neutron contributions. Contributions to the 
transition rates between yrast states do on the other hand mainly stem from the neutron degrees of freedom. In 
conclusion, it seems that the Oj states in ^^Zr and ^"'Zr and the Oj^ state in ^^Zr contain strong collective components 
not reproduced by the present shell model calculation. 
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B. Odd isotopes 



It is somehow surprising to notice that the shell model gives a much better description of the odd than of the 
even Zr isotopes. The reproduction of the low-lying positive parity states are overall satisfactory. On the other hand 
the shell model has some problems in describing the negative parity states. Several of the negative parity states in 
91,932j. g^j-g calculated up to 1 MeV too low. Only a few negative parity states are known in ^^'^^Zr. Thus a detailed 
comparison is difhcult, but the 11/2^ state is reproduced in nice agreement with experiment. 

We will make a detailed study of ^'^Zr. This nucleus is not too simple and not too complex (two protons and three 
neutrons outside the closed core) , and useful information can therefore be extracted from a few central and relatively 
simple configurations. 

The shell- model calculation provides three states below 600 keV, 5/2^'", and 9/2^. From experiments, only 

two states are known, 5/2^ and 3/2^. There are however theoretical arguments supporting a low-lying 9/2^-state. 
The configuration that requires the least energy has all particles in the lowest possible single-particle orbital. In the 
case of ^^Zr the two protons occupy the lpi/2 single-particle orbital, and couple to J = 0. The three neutrons occupy 
the l(i5/2-orbital {\d^i^)^. The three neutrons can then couple to = 3/2+, 5/2+ or 9/2+. Such states will be 

located well below 1 MeV. The lowest experimental 9/2+ candidate observed up to now is seen at 1.46 MeV. 

As many as four 1/2+ states are observed within a small energy interval of 250 keV in the region from 0.95 MeV 
to 1.22 MeV. This observation has no shell-model counterpart. Our calculation provides only one 1/2+ level at 1.40 
MeV. 

We already pointed out that our model has difficulty in describing the negative parity states. Consider for example 
the structure of the ll/2j~ state in ^^Zr, which comes 500 keV lower than the experimental position. Odd parity states 
are constructed by configurations with an odd number of particles in negative parity states. Within our model, protons 
can occupy the lpi/2-orbital and neutrons can occupy the O/iii/2-orbital to produce negative parity states. Both the 
first and second excited 11/2^ states in ^'^Zr are predominantly based on the proton configuration (1^1/2059/2)77- The 
same is true for the ll/2j~ state in the other Zr isotopes. 

Finally, we examine ^^Zr. As we would expect for a system caught in between two "magic" nuclei we recognize 
pronounced single-particle structure. From the occupation numbers, listed in Table VI, we see that both the ground 
state, 1/2+, and the next level, 3/2+, are pure one-quasiparticle states built on a full ld5/2-orbital, i.e. ^^Zr-core. 
Also the 1/2'^ state is a one-quasiparticle state with the l(i5/2-orbital nearly closed, though its calculated position is 
about 0.7 MeV too high. This means that a somewhat lower single-particle energy eogg^^might be more appropriate 
for our effective interaction. The 5/2]^ state can be regarded as a l(i5/2-hole state relative to the ^^Zr-core. All in 
all the yrast states apart from 7/2+ are very well reproduced. In the other, non-yrast states, the ^^Zr-core breaks up 
and two neutrons are distributed equally among the 2sx/2 and 1^3/2 orbitals. 



C. Binding energies 

The binding energies calculated by the formula 

BE(90+"Zr) = BE(^°+"Zr) - BE(««Sr) 

- n(BE(«9Sr) - BE(««Sr)) 

- 2(BE(*^9Y) - BE(8«Sr)) (10) 

are plotted in Figure | In Eq. (|l^) n is the number of valence neutrons. The experimental binding energies show a 
parabola structure with a minimum at ^^Zr, whereas the calculated binding energies increase linearly down to ^®Zr. 
With increasing neutron number the systems become far too strongly bound. This phenomenon of overbinding of 
nuclear systems when effective interactions from meson theory are used has been much discussed in the literature, 
for example in Ref. | pl[ . The solutions to the problem has been that such matrix elements must be modified in order 
to reproduce the binding energies correctly. The so-called centroid matrix elements should be corrected in order to 
reproduce experiment. However, there is no well-defined recipe for doing this. 

The curve labeled "no pn-int" in Figure ^ shows the binding energies with the proton-neutron interaction switched 
off. Now, the systems are too weakly bound, which tells us that the proton-neutron part of the interaction contributes 
strongly to the overbinding. We have made the pn-interaction less attractive by adding an overall constant to the 
diagonal matrix elements. This will not affect the excitation energies relative to the ground state. The constant 
is chosen such as to fit the experimental binding energy of ^°Y. Thus, a constant 0.3 MeV added to the original 
diagonal proton- neutron matrix elements, K^'bablP'^) ~ ^abab(P'^) + 0.3MeV, gives binding energies as shown in Figure 



6 



^, labeled "modified pn-int.". The fit to the experimental values are much improved, although there is a linear rather 
than parabolic dependence on the particle number. 



IV. CONCLUSIONS 



In this work we have performed a full (1^1/20^9/2) proton and (2sld0g7/20/iii/2) neutron shell- model calculation 
of the zirconium isotopes ranging from N=52 to N=60. For the first time we present results from calculations with a 
proton-neutron effective interaction in such heavy nuclei. 

We have succeeded in obtaining a qualitative reproduction of important properties, although there are also short- 
comings. The odd isotopes are very well described by our shell model, in fact better than the calculated even isotopes. 
Both for the odd and the even nuclei we have difficulties in reproducing the negative parity states well. 

For comparison we have presented shell- model results for the neighbouring strontium isotopes. The qualitative 
features are much better reproduced in Zr than in Sr. For example, the shell model fails in reproducing the very 
stable 0'^ — 2+ spacing in Sr. Differently from Zr, there is no sign of N=56 and N=58 sub-shell closures in Sr. The 
quality of the closed-shell core (^^Sr) may in fact be different for the two sets of isotopes. It is likely that the additional 
protons in Zr give a more balanced system and serve to stabilize the core. In Sr the core seems to be softer and more 
unstable. 

The empirical Zr spectra can to a certain extent be interpreted in a weak coupling scheme. The isotopes ^^Zr and 
^^Zr are well described in terms of this model, which in turn indicates that the proton-neutron interaction should not 
be too strong. The simple weak coupling picture collapses however in ^^Zr. 

In order to obtain results for ^^Zr, we performed calculations that are extremely heavy and time consuming. All 
the efforts gave final results that were far off, and in fact a much simpler calculation based on a ^''Sr-core provided 
results in better agreement with the experimental data. 

The occupation numbers give a hint that the 037/2 and the O/111/2 neutron orbitals are of minor importance. The 
major properties are in fact fairly well described within a reduced basis (lpi/2, 0(/g/2)7r (2si/2, 1^5/2, ^d^/2)v- 

In order to further test the wave functions we calculated E2 transition rates in the even Zr isotopes. Transitions 
between yrast states are fairly well reproduced, whereas transitions involving certain excited 0"*" states are calculated 
far too small, indicating that these states contain strong components not accounted for by the present shell model. 

As in other mass regions we fail in reproducing bulk properties such as the binding energies. With increasing 
number of valence nucleons the systems become far too strongly bound. We have demonstrated that this problem 
can be cured by simple adjustments of some selected matrix elements. 

The calculations have been carried out at the IBM cluster at the University of Oslo. Support for this from the 
Research Council of Norway (NFR) is acknowledged. 



APPENDIX A: TABLES 

The low-lying experimental and calculated energy-levels in ^o-s^Zr are listed in Tables g and p| 



[1] F. Schussler, J. A. Pinston, E. Monnand, A. Moussa, G. Jung, E. Koglin, B. Pfeiffer, R.V.F. Janssens and J. van Klinken, 

Nucl. Phys. A 339, 415 (1980). 
[2] K. Becker, G. Jung, K.-H. Kobras, H. Wollnik and B. Pfeiffer, Z. Phys. A 319, 193 (1984). 

[3] H. Mach, M. Moszynski, R.L. Gill, J.A. Winger, John C. Hill, G. Molhar and K. Sistemich, Phys. Lett. B 230, 21 (1989). 
[4] G. Lhersonneau et ai, Phys. Rev. C 56,1445 (1997). 

[5] G. Lhersonneau et at, Phys. Rev. C49, 1379 (1994); G. Lhersonneau et al, Phys. Rev. C 54,1117 (1996). 
[6] D.H. Gloeckner, Phys. Lett. B 42, 381 (1972); D.H. Gloeckner, Nucl. Phys. A 253, 301 (1975). 
[7] P. Halse, J. Phys. G 19,1859 (1993). 

[8] S.S. Ipson, K.C. McLean, W. Booth, J.G.B. Haigh and R.N. Glover, Nucl. Phys. A 253, 189 (1975). 

[9] N. Auerbach and I. Talmi, Nucl. Phys. 64, 458 (1965). 
[10] T. Engeland, M. Hjorth- Jensen, A. Holt and E. Osnes, Phys. Rev. C 48, 535 (1993). 
[11] A. Holt, T. Engeland, M. Hjorth- Jensen, E. Osnes and J. Suhonen, Nucl. Phys. A 618, 107 (1997). 
[12] J. Suhonen, J. Toivanen, A. Holt, T. Engeland, M. Hjorth- Jensen and E. Osnes, Nucl. Phys. A 628, 41 (1998). 



7 



[13] A. Holt, T. Engeland, M. Hjorth- Jensen and E. Osncs, Nucl. Phys. A 634, 41 (1998). 
[14] M. Hjorth- Jensen, T.T.S. Kuo and E. Osnes, Phys. Reports 261, 125 (1995). 
[15] R. Machleidt, F. Sammarruca and Y. Song, Phys. Rev. C 53, 1483 (1996) 

[16] V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen and J.J. de Swart, Phys. Rev. C 48, 792 (1993). 
[17] R.B. Wiringa, V.G.J. Stoks and R. Schiavilla, Phys. Rev. C 51, 38 (1995). 

[18] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[19] T.T.S. Kuo and E. Osnes, Folded-Diagram Theory of the Effective Interaction in Atomic Nuclei, Springer Lecture Notes 

in Physics, (Springer, Berlin, 1990) Vol. 364. 
[20] R.R. Whitehead, A. Watt, B.J. Cole and I. Morrison, Adv. Nucl. Phys. 9, 123 (1977). 
[21] T. Engeland, M. Hjorth-Jensen, A. Holt and E. Osnes, Phys. Scr. T 56, 58 (1995). 

[22] A. Saganek, V. Meyer, S. Mirowski, M. Oteski, M. Siemeaski, E. Wesolowski and Z. Wilhelmi, J. Phys. G 10, 549 (1984). 
[23] T.P. Cleary, Nucl. Phys. A 301, 317 (1978). 

[24] S. Ramavataram, B. Goulard and J. Bergeron, Nucl. Phys. A 207, 140 (1973). 

[25] P. Hoffmann-Pinther and J.L. Adams, Nucl. Phys. A 229, 365 (1974). 

[26] D.S. Chuu, M.M. King Yen, Y. Shan and T. Hsien, Nucl. Phys. A 321, 415 (1979). 

[27] A. Saha, G.D. Jones, L.W. Put and R.H. Sicmsscn, Phys. Lett. B 82, 208 (1979). 

[28] B.A. Brown, P.M.S. Lesser and D.B. Fossan, Phys. Rev. C 13 (1976) 1900. 

[29] M. Hjorth-Jensen, unpublished. 

[30] L. Schellenberg, B. Robcrt-Tissot, K. Kaeser, L.A. Schaller, H. Schneuwly, G. Pricke, S. Gluckert, G. Mallot and E.B. 

Shere, Nucl. Phys. A 333 (1980) 333. 
[31] A.P. Zuker, Nucl. Phys. A 576, 65 (1994); G. Martinez-Pinedo, A.P. Zuker, A. Poves and E. Caurier, Phys. Rev. C 55, 

187 (1997). 

FIG. 1. Experimental low energy level schemes for the Zr (a) and Sr (b) isotopes from N=50 to N=62. 



FIG. 2. Demonstration of the weak coupling scheme. The Zr energy levels are represented by dashed lines. Experimental 
values are used for all levels energies. 



FIG. 3. Selected energy levels in Zr and Z. Numbers in parentheses are the binding energies of the valence nucleons 
outside a **Sr-core. 



FIG. 4. Selected energy levels in Zr and Z. Numbers in parentheses are the binding energies of the valence nucleons 
outside a ®^Sr-core. 



FIG. 5. Calculated energy levels for the Sr isotopes from N=52 to N=60. The experimental levels are shown in Fig. lb. 



FIG. 6. The 2^ excitation energy, experimental (solid) and calculated, complete (dashed) and without pp- and pn-interaction 
(dotted), respectively. 



FIG. 7. Theoretical and experimental spectra of Zr. 



FIG. 8. Relative binding energies ofSOZr-'^^Zr. The curves show the experimental and theoretical binding energy curves, the 
full shell-model calculation, calculations without the pn-interaction and calculations with modified pn-interaction, respectively. 
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TABLE I. Number of shell-model basis states in m-scheme SD representation. 



Tip ~\~ Thji 


Dimension 


Tip ~\~ Thji 


Dimension 


Tip ~\~ Thji 


Dimension 


2 + 


8 


2 + 3 


15 868 


2 + 6 


2 428 814 


2 + 1 


186 


2 + 4 


107 060 


2 + 7 


8 648 777 


2 + 2 


1 572 


2 + 5 


564 393 


2 + 8 


26 201 838 



TABLE IL Single-particle energies, all entries in MeV. 



Single-proton enerf 


;ies 






Single-neutron energies 






jn lpi/2 






1^5/2 


2Sl/2 1^3/2 


0^7/2 


O/111/2 


e{jn) 0.00 


0.90 




0.00 


1.26 2.23 


2.63 


3.50 







TABLE in. 


Occupation 


numbers of the 2f , 4f , 3 


and 5j 


states in ^^-^^Zr. 








O59/2 


lPl/2 


O/111/2 


O57/2 


1^5/2 


1^3/2 


2Sl/2 


«^Zr: 






















0.27 


1.73 


0.02 


0.01 


1.92 


0.02 


0.04 






0.24 


1.76 


0.01 


0.01 


1.96 


0.02 


0.01 




3r 


0.98 


1.02 


0.05 


0.03 


1.67 


0.08 


0.18 




5r 


1.00 


1.00 


0.05 


0.04 


1.79 


0.07 


0.05 


«4Zr: 




















2^ 


0.15 


1.85 


0.06 


0.01 


3.79 


0.06 


0.08 






0.13 


1.87 


0.05 


0.01 


3.87 


0.05 


0.02 




3r 


0.99 


1.01 


0.08 


0.05 


3.52 


0.15 


0.20 




5r 


1.00 


1.00 


0.11 


0.08 


3.48 


0.18 


0.14 


^«Zr: 


















2^ 


0.12 


1.88 


0.09 


0.03 


4.75 


0.13 


1.01 






0.14 


1.86 


0.09 


0.04 


4.78 


1.02 


0.08 




3r 


0.97 


1.03 


0.14 


0.76 


4.52 


0.20 


0.39 




5r 


0.98 


1.02 


0.16 


0.09 


5.21 


0.25 


0.29 


8«Zr: 




















2^ 


0.10 


1.90 


0.12 


0.05 


5.76 


1.10 


0.97 




4^ 


0.18 


1.82 


0.15 


0.98 


5.66 


0.21 


1.01 




3r 


0.99 


1.01 


0.15 


0.55 


5.55 


1.13 


0.63 




5r 


0.91 


1.09 


0.28 


0.21 


5.58 


0.36 


1.57 
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TABLE IV. Occupation numbers of the three lowest-lying 0"*" states in ^®Zr. 





059/2 


lPl/2 


O/111/2 


057/2 


1(^5/2 


1(^3^2 


2si /9 


ZjY. 






















0.33 


1.67 


0.06 


0.03 


1.81 


0.07 


0.03 




0+ 
'-'2 


1.70 


0.30 


0.05 


0.15 


1.61 


0.11 


0.08 




0+ 


1.30 


0.70 


0.07 


0.14 


0.63 


0.06 


1.11 


947-.. 

ZjY. 






















0.17 


1.83 


0.10 


0.04 


3.69 


0.11 


0.06 




0+ 


1.78 


0.23 


0.13 


0.49 


2.75 


0.31 


0.33 




0+ 
^3 


0.42 


1.58 


0.08 


0.06 


2.09 


0.12 


1.65 


9«Zr: 




















0^ 


0.09 


1.91 


0.12 


0.03 


5.66 


0.12 


0.07 




ot 


0.22 


1.78 


0.13 


0.09 


3.88 


0.15 


1.76 




ot 


1.76 


0.24 


0.19 


0.96 


3.86 


0.47 


0.52 


«8Zr: 




















ot 


0.09 


1.91 


0.15 


0.06 


5.76 


0.17 


1.87 






0.44 


1.56 


0.27 


0.72 


5.50 


1.26 


0.26 




03+ 


1.41 


0.59 


0.19 


1.46 


4.79 


0.98 


0.58 



TABLE V. E2 transition rates. Numbers in parentheses indicate uncertainties in the last digit of the quoted experimental 
values. The proton and neutron effective charges are set equal to zero in column 5 and 6, respectively. All entries in Weisskopf 
units (W.U.). 





TRANSITION 


EXP. 




CALC. 






B{E2 








e;" = 1.8,er = 1.5 


= 0.0 


er = 0.0 


'^"Zr 


B{E2 




^0+) 


5.37 (20) 


1.88 








B{E2 




^0+) 


5.2 (10) 


6.21 








bIe2 




^6+) 


2.40 (16) 


2.57 






«^Zr 


B{E2 




^0+) 


6.4 (6) 


2.97 


2.54 


0.02 




B{E2 




^2+) 


14.3 (5) 


0.30 


0.11 


0.77 




B{E2 




^2+) 


4.04 (12) 


0.40 


0.30 


0.01 




B{E2 




^4+) 


> 0.00098 


0.24 


0.002 


0.20 


i^^Zr 


B(E2 


2^ 


^0+) 


4.4 (5) 


3.69 


3.56 


0.001 




B{E2 


ot 


^2+) 


9.3 (4) 


0.01 


0.05 


0.11 




B{E2 




^2+) 


0.876 (23) 


1.26 


1.22 


0.0004 


««Zr 


B{E2 




^0+) 


4(3) 


0.05 


0.03 


0.002 




B{E2 




^0+) 


> 0.020 


0.13 


0.10 


0.002 




B{E2 




^0+) 


> 2.7 


1.17 


1.03 


0.005 


'^^Zr 


B{E2 


2^ 




> 0.24 


0.02 


0.01 


0.0005 




B{E2 


2^ 




> 0.04 


0.64 


0.48 


0.01 




B{E2 




^2+) 


51 (5) 


«0.00 


0.02 


0.02 
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TABLE VI. Occupation numbers in ^''Zr. 





0^9/2 


lPl/2 


O/111/2 


0^7/2 


1(^5/2 


1^3/2 


2Sl /2 


^'Zr: 1/2+ 


0.09 


1.91 


0.11 


0.03 


5.73 


0.15 


0.99 


3/2+ 


0.10 


1.90 


0.11 


0.04 


5.72 


1.04 


0.08 


5/2+ 


0.11 


1.89 


0.12 


0.05 


4.83 


0.14 


1.87 


7/2+ 


0.26 


1.74 


0.15 


1.00 


5.51 


0.17 


0.17 


5/2+ 


0.15 


1.85 


0.10 


0.05 


4.77 


1.06 


1.04 


V2t 


0.14 


1.86 


0.11 


0.05 


4.82 


1.04 


0.99 


3/2+ 


0.18 


1.82 


0.11 


0.06 


4.77 


1.08 


0.98 


ll/2i 


0.88 


1.12 


0.26 


0.07 


5.83 


0.27 


0.88 



TABLE VIL Theoretical and experimental energy levels in Zr and Zr. Energies are given in MeV. 



"Zr 8321. 





SM 




EXP 


J" 


SM 




EXP 


5/2+ 


0.0 


5/2+ 


0.0 


5/2+ 


0.0 


5/2+ 


0.0 


1/2+ 


1.477 


1/2+ 


1.205 


3/2+ 


0.182 


3/2+ 


0.269 


5/2+ 


1.717 


5/2+ 


1.466 


9/2+ 


0.568 


1/2+ 


0.947 


7/2+ 


2.015 


7/2+ 


1.882 


1/2+ 


1.307 


1/2+ 


1.018 


9/2+ 


2.094 


3/2+ 


2.042 


7/2+ 


1.686 


1/2+ 


1.169 


3/2+ 


2.361 


(9/2) + 


2.131 


3/2+ 


1.755 


1/2+ 


1.222 


7/2+ 


2.383 


7/2+ 


2.201 


5/2+ 


1.895 


3/2+, 5/2+ 


1.425 


13/2+ 


2.403 


(5/2,7/2) 


2.367 


5/2+ 


2.114 


(l/2+,3/2,5/2+) 


1.450 


5/2+ 


2.405 


3/2+, 5/2+ 


2.535 


3/2+ 


2.208 


9/2+, 7/2+ 


1.463 


11/2+ 


2.444 


1/2+ 


2.558 


7/2+ 


2.211 


(1/2+, 3/2,5/2+) 


1.470 


5/2- 


1.063 


(11/2)- 


2.170 


13/2- 


1.417 


(9/2-,ll/2-) 


2.025 


15/2" 


1.173 


(5/2)- 


2.190 


9/2- 


1.445 


(9/2-,ll/2-) 


2.363 


13/2" 


1.423 


(13/2)- 


2.260 


11/2- 


1.507 


(9/2-,ll/2-) 


2.662 


11/2" 


1.430 


(15/2)- 


2.288 










7/2" 


1.520 


(11/2)- 


2.321 













TABLE VIIL 


Theoretical and experimental 


energy levels in ^^Zr and ^^Zr 


Energies are given in MeV. 








s^Zr 








s^Zr 




J" 


SM 


J" 


EXP 




SM 


J'^ 


EXP 


5/2+ 


0.0 


5/2+ 


0.0 


1/2+ 


0.0 


1/2+ 


0.0 


1/2+ 


1.021 


1/2+ 


0.954 


3/2+ 


1.062 


3/2+ 


1.103 


3/2+ 


1.285 


3/2+.. 5/2+ 


1.14 


5/2+ 


1.168 


7/2+ 


1.264 


7/2+ 


1.482 


3/2+, 5/2+ 


1.324 


7/2+ 


1.940 


(5/2+) 


1.400 


5/2+ 


1.613 


7/2+, 9/2+ 


1.618 


5/2+ 


2.427 


(5/2+) 


1.859 


9/2+ 


1.857 


(3/2)+ 


1.618 


7/2+ 


2.457 


(5/2+) 


1.997 


3/2+ 


2.051 


(5/2)+ 


1.722 


3/2+ 


2.604 


(3/2,5/2) 


2.058 


7/2+ 


2.306 


1/2^+^3/2,5/2+ 


1.904 


9/2+ 


2.644 


(7/2)+ 


2.234 


1/2+ 


2.414 


l/2(+),3/2,5/2+ 


1.940 


5/2+ 


2.759 


(7/2)+ 


2.508 


5/2+ 


2.490 


5/2(+) 


1.956 


7/2+ 


2.786 


(7/2)+ 


3.161 


13/2- 


1.919 


9/2-, 11/2- 


2.025 


11/2- 


2.356 


(7/2-) 


1.807 


11/2- 


2.039 


l/2-,3/2- 


2.816 


9/2- 


2.492 


(11/2-) 


2.264 
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TABLE IX. Low-lying states in ^°Zr and ^^Zr. Energies are given in MeV. 



8°Zr ''^Zr 



r 


**Sr-core 


EXP 




^°Zr-core 


*®Sr-core 


EXP 


ot 

2 


1.706 


1.761 


1 


2.738 


1.693 


1.383 


2t 


2.003 


2.186 




0.601 


0.581 


0.935 


2+ 




3.309 


2+ 


1.792 


1.921 


1.847 








2+ 


3.093 


2.167 


2.067 


3r 




2.748 


3r 


4.018 


1.904 


2.340 




2.235 


3.077 


4^ 


0.850 


0.823 


1.496 


4r 


1.693 


2.739 


42+ 


2.677 


2.379 


2.398 






2.319 


5r 


4.261 


1.316 


2.486 


6^ 


2.317 


3.448 


6^ 


3.206 


2.596 





TABLE X. Low-lying states in ^*Zr and ^^Zr. Energies are given in MeV. 



*Zr 



"Zr 



Zr-core 


®®Sr-core 


EXP 


r 


90 rv 

Zr-core 


*'^Sr-core 


EXP 


2.254 


2.213 


1.300 


02+ 


1.602 


2.228 


1.582 


0.557 


0.520 


0.919 


2^ 


1.097 


1.426 


1.751 


1.490 


1.764 


1.671 


2^ 


2.105 


2.661 


2.226 


1.930 


2.228 


2.151 


2+ 


2.211 


2.768 


2.669 


3.923 


2.223 


2.058 


3r 


3.833 


3.732 


1.897 


0.889 


0.817 


1.470 


4f 


2.054 


2.528 


2.750 


1.898 


2.162 


2.330 




2.230 


2.869 


2.857 


4.200 


2.673 


2.945 


5r 


4.141 


2.306 


3.120 



3r 

4^ 

5r 





TABLE XI. Low-lying 


states in Zr. Energies are 


given in MeV. 






'"'Zr-core 


»«Zr 
^^Sr-(:(jrc 


''\Sr-core 


EXP 


02+ 

0^ 
2^ 
2J 
2+ 

3r 

4^ 
4J 

5r 


1.641 
2.773 
1.300 
2.207 
2.442 
3.983 
2.147 
2.048 
3.576 


1.904 
2.471 
1.463 
2.619 

3.579 
2.449 
2.650 
3.478 


0.529 
1.969 
1.152 
1.773 
2.149 
2.597 
1.574 
1.918 
1.318 


0.854 
1.859 
1.223 
1.591 
1.744 
1.806 
1.843 
2.330 
2.800 
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> 2 



1- 



5- 

2+ 



0+ 



0+ 



90 



Zr 



2+ 
4+ 



92 



Zr 



2+ 



2+ 



2+ 



94 



Zr 



96 



Zr 



3- 

2+ 
2+ 
0+ 
.2+ 



^Zr 



(a) 



6+ 
2+ 

4+ 
0+ 



100 



Zr 



6+ 
(0+) 

4+ 

2+ 
0+ 



102 



Zr 



3- 



2+ 
0+ 



2+ 



1- 



0- 0+ 



^Sr 



2+ 
4+ 



90 



Sr 



■ 2(+). 

■ 2+ 



.(0+) 
■ 3- 

-(0+) 



92 



Sr 



94 



Sr 



(2+) 
2+ 
0+ 



96 



Sr 



(b) 



.(2+) 
(6+) 

■ 4+ 
. 0+ 



^Sr 



8+ 



(6+) 



4+ 

2+ 
0+ 



100 



Sr 



1 



A + 



2,3 

3~ 
2+ 



0+ 
4+ 



. 4+ 

■ 3" 
■(2 + ) 

■ 2 + 

■ 2 + 

■ 4+ 

■ 0+ 



A + 



■ 2i 

■ 4+ 
■(2+) 

■ 3" 



2 + 
4+ 
0+ 



A + 



+ 



■ (2+,3+,4+) 

3<"' 

■ 0+ 

(0+) 



■ 0+ 

■ 4± 

■■ 0+ 

(2 + ) 

■ 3+ 

■ 2+ 



3~ 
2+ 
0+ 



2+ 



(2+) 



^Zr 



^Sr + «"Zr 



*Zr 



'Sr + «"Zr 



"Zr 



1 



2+ 

(-4.39) (-4.30) (-6.62) (-8.23) 

0+ 0+ 0+ 0+ 

SM EXP SM SM EXP SM 

''°Zr-core ^^Sr-core ^°Zr-core ^^Sr-core 



1 



5- — . 962r ^ 3- . 982j 



3- 



5- 



/ ■. ■• 5- ; / 4+ 

4+ ■■ / 0+ 4+ : :/ 

/ / 0+ 

2+ ^ 



\ 

\ / 



/ 



(-8.21) (-12.36) (-7.48) (-15.16) 

0+ 0+ 0+ 0+ 

SM EXP SM SM EXP SM 

''°Zr-core ^^Sr-core ^°Zr-core ^^Sr-core 



1 




1 




1 




1 




1 



